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1.  Introduction 


Brain  tumors  are  most  frequent  solid  cancer  among  all  kinds  of  childhood  cancer. 
Heterogeneity,  invasiveness,  and  complex  microenvironment  are  making  therapies  as 
well  as  research  on  asrtocytic  brain  tumor  particularly  challenging  Low-grade  gliomas 
are  curable  and  most  patients  could  live  without  further  progression  and  severe  condition 
for  many  years.  However,  once  the  glioma  progress  to  high  grade,  the  quality  of  life  and 
survival  of  patients  are  very  poor  '.  Based  on  previous  work,  we  hypothesized  the  bone 
marraw  derived  cells  (BMDCs)  could  serve  as  a  mediator  of  transition  of  low  grade  to 
high-grade  tumor  The  proposed  project  aims  to  characterize  the  various  lineages  of 
glioma  associated  BMDCs  including  hematopoietic,  endothelial,  and  mesenchymal 
lineages  in  both  low  grade  and  high  grade  stages  of  glioma.  Studying  the  effect  on  glioma 
transition  from  low  grade  to  high  grade  by  depleting  distinct  populations  of  bone  marrow 
derived  inflammatory  cells  including  monocyte,  granulocytes,  endothelial  progenitors, 
and  mesenchymal  progenitors.  Dissecting  molecular  mechanism/signaling  of 
differentiation  of  glioma  associated  BMDCs,  and  screening  the  key  factors  or  targets 
through  the  entire  regulatory  pathway.  It  would  contribute  to  develop  therapeutic 
strategies  to  target  a  specific  population  of  BMDCs  and  their  subsequent  recruitment,  in 
order  to  suppress  the  malignant  transformation  of  gliomas.  In  this  project,  we  have 
initiated  the  study  of  BMDCs  with  RCAS  and  G1261  murine  glioma  models  as  well  as 
glioma  patients.  We  also  used  transgenetic  tools  to  deplete  certain  populations  of  BMDCs 
to  study  functional  contribution  of  BMDCs  for  glioma  progression. 
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2.  Keywords 
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3.  Overall  Project  Summary 


The  project  focuses  on  studying  the  microenvironment  and  functions  of  bone  marrow 
derived  cells  within  pediatric  astrocytic  tumor.  In  this  study,  two  major  glioma  models 
will  be  used  to  investigate  the  role  of  BMDCs  primarily.  One  is  a  transgenic  mouse  using 
the  RCAS/Tv-a  system  created  by  Holland  and  Varmus  that  develops  low-grade  gliomas 
which  progress  to  high-grade  tumors  over  the  course  of  twelve  weeks  The  other 
model  is  syngeneic  orthtopic  glioma  model.  G1261  a  C57/B16  derived  glioma  cell  line 
with  different  markers  such  as  GFP  or  luciferase,  was  intracranial  injected  to  C57/B16 
mice  to  create  allograft  glioma.  In  addition,  xenograft  models  with  human  glioma  cell 
lines  are  also  utilized.  Furthermore,  we  also  have  used  glioma  patients’  blood  samples  to 
analyze  various  lineages  of  BMDCs 

Task  1.  Charactering  the  various  lineages  of  glioma  associated  BMDCs  including 
hematopoietic,  endothelial,  and  mesenchymal  lineages  in  both  low  grade  and  high 
gradestages  of  glioma.  (90%  complete) 

We  have  analyzed  the  myeloid  lineage  and  endothelial  lineage  of  BMDCs  in  both 
patients’  blood  samples  and  murine  glioma  models  at  low-grade  or  high-grade  stage. 

la.  Assess  the  frequency  and  absolute  numbers  of  circulating  BMDCs  in  glioma 
patients.  Using  flow-cytometry  of  hematopoietic,  as  well  as  endothelial  and 
mesenchymal  markers,  we  will  investigate  whether  the  frequency  of  HPCs, 
VEGFR2+  EPCs  and  CD105+  MSCs  correlates  with  glioma  transformation  in 
patients.  Human  subjects  involve  this  study  will  be  radiographically-suspected  or 
biopsy  prove  low-grade  astrocytoma  (WHO  grade  I  and  H)  versus  histology  proven 
high  grade  gliomas  (WHO  grade  III  and  IV).  Blood  samples  will  be  collected  at  the 
time  of  diagnosis.  Subjects  selected  for  this  study  may  be  between  2  years  and  16 
years  of  age.  We  plan  to  enroll  160  human  subjects,  and  quarterly  enrollment  is  20. 
The  identifications  of  human  subjects  will  not  be  accessible  for  research  team,  and 
only  information  of  human  subjects  that  research  team  is  aware  of  is  the  patients’ 
diagnosis  and  diseases’  history. 

Up  to  date,  we  have  recruited  107  glioma  patients  with  low-grade  or  high-grade  glioma, 
plus  26  healthy  volunteers  as  control  into  this  study.  Patients  or  healthy  volunteers’ 
peripheral  blood  has  been  analyzed  with  various  lineages  makers  as  showed  in  the  figures 
below: 
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Figure  1.  Characterizing  myeloid  lineage  of  BMDCs  in  patients  (CDllb,  CD33,  CD14, 
and  CD  16)  by  flow  cytometery  in  peripheral  of  low-grade  astrocytoma  patients  (LGA)  vs 
glioblastoma  patients  (GBM). 
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Figure  2.  IHC  of  CDllb  (infiltrated  myeloid  cells)  on  archived  paraffin  embedded  tumor 
tissue  from  low-grade  astrocytoma  patients  (grade  II)  vs  glioblastoma  patients  (grade  IV). 
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Figure  3.  Characterizing  endothelial/myeloid  lineage  of  BMDCs  in  patients  by  CD  11b, 
KDR  (VEGFR2)  in  peripheral  of  low-grade  astrocytoma  patients  (LG A),  glioblastoma 
patients  (GBM),  and  healthy  volunteers. 


Figure  4.  Statistic  dots-plot  on  frequency  of  CDllb-i-/VEGFR2(KDR)-i-  cells  out  of  total 
CDllb-i-  cells  in  patients  with  different  stages  of  disease.  Healthy  donor  served  as  a 
control.  GIII/GIV  vs  healthy,  One  way  ANOVA,  P  <0.001. 
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Eigure  5.  Grade  II  glioma  patients  were  sub-divided  into  two  groups,  based  on  the 
frequency  of  CDllb-t/KDR  (VEGER2)-l-  cells:  KDR  Hi  and  KDR  Lo  (cut-off,  1.0%). 
Kaplan-Meier  curves  for  progression-free  survival  over  12  months  are  presented  (27 
subjects).  P  <0.01  by  log-rank  test. 

During  the  study  of  BMDCs  with  glioma  patients’  blood  samples,  we  have  found  that  the 
number  of  myeloid  derived  suppor  cells  (MDSCs)  within  myeloid  lineage  increased 
following  the  progression  of  diseases.  The  MDSCs  are  heterogeneous  regarding  the 
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expression  of  CD  14/CD  16,  representing  monocytic  or  granulocytic  sub-lineages  (Figure 
1).  While  the  tumor  progressed,  we  also  observed  more  infiltrated  myeloid  cells  within 
tumor  tissues  (Figure  2).  Interestingly,  when  we  examine  the  expression  of  VEGFR2  to 
study  endothelial  lineage  of  BMDCs,  we  found  the  majority  of  VEGFR2  expressing  cells 
are  CDl  lb  positive  as  showed  in  Eigure  3.  It  suggested  that  there  is  overlap  or  interaction 
between  myeloid  lineage  and  endothelial  lineage  of  BMDCs.  Quantification  of 
CD105-l-/C-kit-l-  BMDCs  indicates  no  significant  difference  in  MSCs  from  patients  with 
in  low-grade  patients  vs  high-grade  patients.  Therefore,  we  would  focus  on  raveling  the 
role  of  myeloid/endothelial  lineage  in  glioma  progression  by  following  animal  study. 

Based  on  the  analysis  on  large  cohort  of  patients,  we  found  that  the  number  of 
CDllb-t/KDR-l-  BMDCs  have  significant  correlation  with  the  pathological  diagnosis  of 
patients.  The  ratio  of  KDR-l-  myeloid  cells  was  significantly  higher  in  blood  samples 
derived  from  patients  with  high-grade  gliomas  when  compared  with  samples  from 
patients  with  low-grade  gliomas.  Specifically,  the  levels  of  CDllb-i-KDR-i-  cells  were 
higher  in  patients  with  WHO  grade  III  or  IV  tumors  than  it  is  in  patients  with  grade  I 
tumors  and  healthy  controls  (Eigure  4).  Perhaps  most  striking  is  the  apparent  delineation 
of  two  distinct  populations  of  patients,  each  with  pathologically  defined  grade  II 
astrocytoma.  We  sub-divided  patients  with  grade  II  astrocytoma  into  two  subgroups, 
those  with  high  versus  low  CDllb-l-KDR-l-  values  (cut-off  of  1.0%).  The  clinical  course 
and  medical  histories  of  the  patients  were  closely  followed  over  12  months.  Disease 
progression  in  patients  with  grade  II  glioma  was  assessed  by  MRI  and  through 
subsequent  histological  diagnoses  of  grade  III  or  IV  glioma.  The  progression-free 
survival  of  our  cohort  of  grade  II  patients  is  shown  in  Eig.  2c.  Patients  with  an  elevated 
percentage  of  CDllb-l-KDR-l-  (KDR  Hi)  myeloid  cells  presented  a  significantly  higher 
likelihood  of  tumor  progression  from  fibrillary  to  anaplastic  variants  (Eigure  5). 


lb.  Evaluate  the  frequency  and  number  of  the  same  BM-derived  populations  in 
blood,  bone  marrow  and  metastatic  organs  of  murine  models  of  glioma  during  low 
grade,  transformation  grade  and  high  grade  phases.  Specifically,  we  will  investigate 
the  mobilization  of  HPCs,  EPCs,  and  MSCs  by  flow  cytometry. 

We  have  analyzed  the  HPCs,  EPCs,  and  MSCs  in  both  RCAS  and  G1261  murine  glioma 
models.  The  4~6  weeks  post-injection  of  RCAS  model  were  consider  as  the  low-grade 
stage,  and  6~9  weeks  were  considered  high-grade  stage.  The  G1261  model  was 
considered  as  high-grade  glioma  model. 
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Week  5  Week  8 

Figure  6.  RCAS-tva  murine  glioma  model  on  low-grade  stage  (week  5)  and  high-grade 
stage  (week  8).  The  MRI  or  H&E  staining  from  representative  mice  were  showed. 


Figure  7.  The  Kaplan-Miere  Survival  curves  of  RCAS-tva  tumor  bearing  mice  with  or 
without  bone  marrow  transplantation  (BMT).  The  median  survival  time  for  RCAS  mice 
without  BMT  is  7  weeks,  and  with  BMT  is  1 1  weeks. 

We  characterized  the  RCA-tva  murine  glioma  model  in  our  experimental  setting,  and 
figured  out  their  low-grade  stage  and  high-grade  stage  evidenced  by  MRI  and  histology. 
Additionally,  we  showed  the  RCAS  mice  with  bone  marrow  transplantation  have  delayed 
progression  of  tumor,  which  is  important  for  our  next  step  of  study. 
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Figure  8.  The  expression  of  VEGFR2  on  BMDCs  in  G1261  model.  Upper  panel  of 
schematic  model  showed  the  VEGFR2-GFP  knock  in  mice  for  studying  expression 
pattern  of  VEGFR2.  Lower  panel  of  flowcytometry  graphs  indicate  the  expression 
VEGFR2  on  GDI  lb+  or  GDI  lb-  population. 


Non-tumor  mice 
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Eigure  9.  EPG  in  G1261  glioma  model.  Elowcytometry  graph  of  lineage  negative  GDI 44-1- 
VEGFR2-I-  EPG  in  G1261  tumor  bearing  mice. 
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We  have  studied  the  myeloid  lineage  and  EPC  with  G1261  glioma  model.  We  observed 
that  both  VEGFR2+  myeloid  cells  and  EPCs  were  elevated  by  G1261  tumor,  indicating 
potential  interaction  of  BMDCs  differentiation  with  primary  glioma. 


Task  2.  Studying  the  effect  on  glioma  transition  from  low  grade  to  high  grade  by 
depleting  distinct  populations  of  bone  marrow  derived  inflammatory  cells  including 
monocyte,  granulocytes,  endothelial  progenitors,  and  mesenchymal  progenitors. 
(90%  complete) 

We  have  set  up  all  trans-genetic  mice  lines  including  ITGAM(CDllb)-DTR/EGFP  mice, 
RosaCreERT2/PDGFRaloxP/loxP  mice,  and  RosaCreERT2WEGFR21oxP/loxP  mice  in 
suitable  genetic  background  for  bone  marrow  transplantation  experiments.  We  performed 
the  lineage  depletion  experiments  with  ITGAM(CDllb)-DTR/EGFP  and 
RosaCreERT2WEGFR21oxP/loxP  mice  for  myeloid  or  endothelial  lineages  in  both 
RCAS  and  G1261  tumor  models. 


2a.  We  plan  to  transplant  the  bone  marrow  from  ITGAM(CDllb)-DTR/EGFP  mice 
into  the  RCAS  and  G1261  glioma  models,  and  use  diphtheria  toxin  to  induce 
depletion  of  myeloid  cells  in  this  RCAS  glioma  model. 
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Figure  10.  Competitive  BMT  showed  that  KDR  knockout  in  bone  marrow  cells  lead  to 
deficiency  of  the  differentiation  of  monocytic  cells.  Lethal  dose  irradiated  C57/bl6  mice 
were  transplanted  with  UBC-GFP/rosa26ERT2-cre/KDRfl/fl  and  rosa26ERT2- 
cre/KDRfl/+  bone  marrow  cells,  and  G1261  tumors  were  implanted  after  bone  marrow 
engraftment.  (a)  Peripheral  white  blood  cells  were  analyzed  on  SSC  and  GFP  by  flow 
cytometry.  The  GFP+  and  GFP-  populations  were  gated  for  further  analysis.  CD  11b  vs 
GFP  (b),  Ly6C  vs  Ly6G  (c). 
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We  have  developed  the  ITGAM(CDllb)-DTR/EGFP  bone  marrow  transplanted  mice 
and  implanted  mice  with  G1261  tumor  after  bone  marrow  was  engrafted.  We  have  tried  to 
deplete  the  GDI  lb  positive  cells  in  tumor  bearing  mice  once  with  5  mice  in  each  group. 
However,  we  only  obtained  approximately  10%  deduction  of  GDI  lb  cells  compared  with 
control  group.  In  this  case,  we  didn’t  see  effect  on  tumor  growth.  However,  as  Figure  10 
showed,  knockout  KDR  (VFGFR2)  led  to  significant  reduction  of  CDllb+  cells,  among 
of  which  Ly6  C+  cells  were  suppressed  dramatically.  In  conclusion,  knockout  VFGFR2 
has  considerably  effect  on  myeloid  lineage,  which  is  inconsistent  co-expression  VFGFR2 
and  GDI  lb  in  high-grade  patients’  BMDCs  samples.  Accordingly,  we  have  used 
VFGFR2  deficient  model  in  the  following  study. 

2b.  Bone  marrow  from  RosaCreERT2A^EGFR21oxP/loxP  mice  will  be  transplanted 
to  both  RCAS  glioma  and  G1261  bearing  mice  to  deplete  endothelial  lineage  of  bone 
marrow  derived  cells  by  knocking  out  the  VEGFR2  gene.  The  total  number  of  mice 
will  be  used  is  30. 


We  have  successfully  knockout  VEGFR2  in  both  RCAS  model  and  G1261  murine  model. 
We  have  studied  the  effect  of  VEGFR2  deficiency  on  glioma  progression. 

As  Figure  11  showed,  we  could  obtain  around  90%  knocking  out  efficiency  with  our 
RosaCreFRT2WFGFR21oxP/loxP  system,  without  affecting  counts  of  blood  cells  (CBC) 
of  mice. 
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Figure  11.  Validation  of  depletion  of  VFGFR2  (KDR)  allele  in  bon  marrow  cells,  (a) 
Amplification  of  KDR-l-  or  KDR-  allele  on  bone  marrow  cells  of  mice  bearing  indicated 
genetic  background.  Tamoxifen  was  applied  to  mice  one  week  before  testing.  Complete 
bloods  count  on  mice  with  WT,  KDR-l-/-,  or  KDR-/-  bone  marrow  (b)  Red  blood  cells,  (c) 
white  blood  cells,  and  (d)  Platelets. 
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When  we  knockout  VEGFR2  from  BMDCs  in  G1261  models,  we  observed  that  tumor 
progression  were  suppressed  as  showed  in  figure  12.  In  the  tumor  tissue,  we  found  much 
less  tumor  associated  myeloid  cells  (CDllb+)  in  the  VEGFR2  KO  group  compared  with 
control  group  (Figure  13).  With  RCAS  system,  we  found  the  similar  phenotype,  after  we 
performed  the  bone  marrow  transplantation  and  induced  VEGFR2  knockout,  we 
observed  that  tumor  progression  were  delayed  and  median  survival  time  (MST)  was 
significantly  elongated  (Figure  14). 


05 


Figure  12.  Knocking-out  VEGFR2  (KDR)  in  BMDCs  suppress  G1261  tumor  growth, 
tumor-associated  myeloid  cells,  and  vasculatures.  Chimeric  C57/bl6  mice  transplanted 
with  rosa26ERT2-cre/KDRfl/fl  bone  marrow  cells  (labeled  as  BM-KDR  KO,  and  BM- 
KDR  control  is  rosa26ERT2-cre/KDRfl/-i-)  were  implanted  with  luciferase  labeled-G1261 
tumors  intracranially.  Tamoxifen  were  applied  at  day  3  post-  implantation.  The  tumors 
were  monitored  by  bioluminescence.  The  quantification  of  bioluminescence  based  tumor 
growth. 
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Figure  13.  Immuno-staining  of  vascular  basement  (Collagen  IV)  and  myeloid  cells 
(CD  11b)  in  the  G1261  tumors  from  each  group.  Scale  bar,  20  pm.  The  number  of  CD  11b 
cells  in  tumor  core  or  tumor  margin  was  quantified. 
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Figure  14.  Knocking-out  VEGFR2  (KDR)  in  BMDCs  suppress  spontaneous  malignant 
transformation  of  RCAS  gliomas,  (a)  Flow  chart  of  experimental  design.  Oncogenes  were 
transduced  to  P0-P2  pups  in  the  beginning.  Mice  received  lethal  dose  irradiation  and 
received  bone  marrow  transplantation  at  week  2,  and  then  tamoxifen  was  applied  to  mice 
at  week  5  to  induce  the  ablation  of  target  gene.  Tumors  were  monitored  by  MRI  over  the 
process,  (b)  Kaplan-Meier  symptom  free  survial  curve  for  RCAS  mice  transplanted  with 
rosa26ERT2-cre/KDRfl/fl  bone  marrow  cells  (BM-KDR  KO),  rosa26ERT2-cre/KDRfl/fl 
bone  marrow  cells  (BM-KDR  Con),  or  without  irradiation/transplantation  (W/O  BMT). 
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Week  12 


□  No  Tumor 

□  Low-Grade  Tumor 


■  High-Grade  Tumor 
£□1  Mortality 


Eigure  15.  At  weeks  6,  week  9,  and  week  12,  tumors  in  each  group  of  mice  were  assessed 
and  graded  based  on  their  MRI  and  further  confirmation  with  histology.  N=  19-26. 

We  compared  brain  tumors  in  RCAS/tva  mice  at  weeks  6,  9  and  12,  utilizing  MRI 
followed  by  histological  examination  (Eigure  15).  Knocking  out  KDR  in  BMDCs 
significantly  delayed  the  malignant  transformation  of  low-grade  glioma.  By  week  9,  60% 
of  RCAS/tva  mice  transplanted  with  control  R26CreERKDRfl/-l-  bone  marrow  developed 
high-grade  tumors,  which  had  typical  malignant  features,  including  robust  gadolinium 
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enhancement  on  MRI,  pseudopalisading  necrosis,  and  microvascular  proliferation.  In 
contrast,  tumors  in  70%  of  RCAS/tva  mice  transplanted  with  R26CreERKDRfl/fl  bone 
marrow  remained  in  the  low-grade  stage  without  signs  of  malignant  transformation  at 
week. 

2c.  Bone  marrow  cells  from  RosaCreERT2/PDGFRaloxP/loxP  mice  will  be 
transplanted  to  each  genetic  and  orthotopic  glioma  generating  mice.  PDGFRa  is 
expressed  on  bone  marrow  derived  mesenchymal  stem  cells  and  PDGF-  PDGFRo 
signal  axis  is  very  critical  for  maintenance  of  mesenchymal  lineage.  By  knocking  out 
PDGFRa  gene,  we  will  study  the  influence  of  defection  of  mesenchymal 
differentiation  on  the  progression  of  low-grade  glioma.  The  total  number  of  mice 
will  be  used  is  30. 

We  have  completed  the  crossing  donor  RosaCreERT2/PDGFRaloxP/loxP  mice,  and  the 
bone  marrow  transplantation  experiments  along  with  tumor  study  are  ongoing.  As  study 
has  showed.  Deletion  of  PDGFRa  gene  in  BMDCs  have  no  significant  influence  on 
development  of  G1261  tumor  in  mice  brain.  Interestingly,  we  observed  that  tumor  did 
growly  slower  and  in  RosaCreERT2/PDGFRaloxP/loxP  mice  compared  with  control 
mice.  It  suggests  that  residential  PDGFRa-l-  cells  play  more  important  role  in  tumor 
progression  than  bone  marrow  derived  PDGFR-l-  cells.  (Figure  16) 


★* 

PDGFRa*/*  PDGFRa^/^ 


Figure  16.  Knocking  out  PDGFRa  in  stromal  cells  inhibits  glioma  progression. 
Luciferase-G1261  cells  were  intracranially  injected  into  RosaCre-ERT2/ 

PDGFRa  fl/fl  mice  and  RosaCre-ERT2/PDGFRa-l-/-l-  mice.  A,  tumor  burden  is 
represented  by  photon  flux  intensity.  P  <  0.05.  B,  precontrast  and  postcontrast  Tl- 
weighted  MRIs  were  performed  on  each  group  at  day  15.  Arrows,  enhancing  tumor. 
Tumor  volume  quantification  is  demonstrated.  P  <  0.01.  C,  vWF  was  stained  in  the  tumor 
periphery  and  within  the  tumor  of  each  group  (left  and  middle  panels).  Collagen  IV  was 
stained  in  each  group  (right).  Scale  bar,  50  pm.  Vessel  diameters  were  quantified.  P  < 
0.05. 
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Task  3.  Dissecting  molecular  mechanism/signaling  of  differentiation  of  glioma 
associated  BMDCs,  and  screening  the  key  factors  or  targets  through  the  entire 
regulatory  pathway.  (60%  complete) 

We  will  proceed  to  study  which  genes  in  certain  lineages  of  BMDCs  could  play  a  critical 
role  in  promoting  the  invasion  of  glioma.  We  will  isolate  specific  population  from 
BMDCs  from  the  tumor  tissue,  blood,  and  bone  marrow.  Then  we  utilize  microarray/next 
generation  seuqencing,  antibody  array,  LC-TOF-TOF  to  detect  gene  and  protein 
expression.  To  determine  the  functional  contribution  of  certain  genes  of  interest,  or 
certain  subpopulation  of  BMDCs,  we  plan  to  utilize  a  series  of  in  vitro/in  vivo 
experiments,  including  basement  invasion  assay  and  knock  specific  genes  bone  marrow 
cells  and  then  transplant  them  into  glioma  producing  mice,  to  elucidate  their  specific 
roles  in  promoting  invasion  of  glioma  cells  and  infiltration  of  glioma.  Total  number  of 
mice  will  be  used  for  this  study  is  30. 

We  have  been  working  on  optimizing  sorting  different  lineages  of  BMDCs  by  FACS  or 
MACS  and  tested  a  few  samples  by  RNA  sequencing.  The  data  suggested  that 
ID2WEGF2  signaling  was  playing  important  role  in  myeloid  differentiation. 

In  order  to  further  delineate  the  signaling  network  driving  myeloid/endothelial  lineage 
differentiation,  we  performed  gene  expression  profiling  of  VEGFR2-expressing 
hematopoietic  progenitor  cells  (Lin-C-kit+)  by  mRNA  sequencing.  Differentially 
expressed  genes  from  VEGER2+  versus  VEGER2-  HPCs  were  clustered  and  arranged  in 
a  heatmap.  Differentially  expressed  genes  were  also  clustered  and  displayed  (Eigure  17). 
Candidate  genes  (P  <0.05,  >1. 5-fold  change)  were  divided  according  to  subsets  with  the 
highest  expression  and  analyzed  for  categories  with  significant  enrichment  (P<  0.05)  of 
categories  in  Gene  Ontology  (GO)  biologic  processes  using  DAVID  tools.  Similar 
categories  were  grouped  accordingly  (Eigure  18).  Inhibitor  of  DNA  binding  proteins  2 
(ID2)  was  identified  as  a  significantly  up-regulated  gene  in  VEGER2-I-  HPCs  a  strong 
candidate  to  be  an  upstream  molecule  mediating  myeloid  endothelial  differentiation. 


Bone  Marrow  Hematopietic  Progenitors 


Color  Key 
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Figure  17.  RNA-sequencing  based  gene  profiles  on  KDR+  (VEGFR2)  peripheral  and 
bone  marrow  cells.  The  heatmap  of  significantly  altered  genes  in  Lin-c-kit+Sca-1-  from 
naive  mice,  Lin-c-kit+Sca-1-  KDR-GFP+  cells  and  Lin-c-kit+Sca-l-KDR-GFP-  cells 
from  bone  marrow  of  tumor  bearing  mice. 


I  transforming  growth  factor,  beta  receptor  II 
CD300  antigen  like  family  member  F 
tescalcin,  similar  to  Tescalcin 
ATPase,  transporting,  lysosomal  VO  subunit  A1 
catenin  (cadherm  associated  protein),  beta  1 
avian  reticuloendotheliosis  viral  (v-rei)  oncogene  related  B 
lysosomal  acid  lipase  A 
coagulation  factor  II  (thrombin)  receptor 
peroxiredoxin  2 

B-cell  leukemia/lymphoma  2  related  protein  Ala 

predicted  gene  8566.  superoxide  dismutase  1.  soluble,  similar  to  Superoxide  dismutase 

peroxiredoxin  1  predicted  gene  7204 

transforming  growth  factor,  beta  1 

similar  to  Leukocyte  immunoglobulin-like  receptor 

hemoglobin  alpha,  adult  chain  2;  hemoglobin  alpha,  adult  chain  1 

predicted  gene  4327;  predicted  gene  8683;  similar  to  40S  ribosomal  protein  S19, 

predicted  gene  5618.  predicted  gene  6476;  predicted  gene  4796; 

Yamaguchi  sarcoma  viral  (v-yes  -1)  oncogene  homolog 
hematopoietic  cell  specific  Lyn  substrate  1 
inhibitor  of  DNA  binding  2 
deoxyribonuclease  11  alpha 


Figure  18.  The  candidate  genes  (P  <0.05,  >1.5-fold  change;  Lin-c-kit+Sca-1- 
KDR(VEGFR2)-GFP+  vs  Lin-c-kit+Sca-1 -KDR-GPP-)  were  divided  according  to  the 
subset  with  the  highest  expression  and  analyzed  for  categories  with  significant 
enrichment  (P<  0.05)  of  categories  in  Gene  Ontology  (GO)  biologic  processes  using 
DAVID  tools.  Similar  categories  were  grouped  accordingly. 
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The  Level  of  ID2  in  Hematopoietic  Cells 
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Figure  19.  The  expression  of  ID2  in  various  lineages  of  hematopoietic  cells.  Means  + 
SEM,  ***p<  0.0001  by  one  way  ANOVA. 


The  Level  of  ID2  In  CD45-CD34+  HPC) 


in  Glioma  Patients 


Figure  20.  The  expression  of  ID2  in  hematopoietic  progenitor  cells  (CD45-CD34+)  from 
patients  with  low-grade  or  high-grade  gliomas.  Means  +  SEM,  P<  0.0001  by  one  way 
ANOVA. 

To  validate  the  up-regulation  of  ID2  in  HPCs,  we  performed  real-time  PCR  on 
HPCs  and  LybC^/LybC^  myeloid  cells  from  tumor-bearing  or  naive  mice.  KDR*^^"^ 
HPCs  from  tumor-bearing  mice  had  higher  ID2  levels  than  KDR*^^'  HPCs  and  HPCs 
from  naive  mice  (Figure  19).  The  expression  of  ID2  was  generally  lower  in  myeloid  cells 
than  in  HPCs.  We  also  examined  the  expression  of  ID2  in  glioma  patients.  HPCs  were 
isolated  from  healthy  controls  and  patients  diagnosed  with  low-grade  and  high-grade 
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gliomas.  HPCs  from  high-grade  glioma  patients  had  significantly  higher  levels  of  ID2 
than  did  HPCs  from  healthy  controls  and  low-grade  glioma  patients  (P<0.0001,  Figure 
20). 
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Figure  21.  Upstream  of  KDR  (a)  the  bioinformatic  analysis  of  promoter  and  5’-UTR 
region  of  KDR  by  biobase.  Highlighted  fragment  of  KDR  promoter  showed  potential 
binding  site  for  E2A.  (b)  Chromatin  immunoprecipitation  of  E2A  with  or  w/o  over¬ 
expression  of  ID2.  Probe  on  KDR  promoter  fragments  were  used.  Means  +  SEM, 
**P<0.01,  ***P<0.0001  by  one  way  ANOVA.  (c)  Immunobloting  of  KDR  and  ID2  on  in 
vitro  cultured  Lineage  negative  bone  marrow  cells  treated  with  or  w/o  TGE-pi  and 
GMCSF.  Phospho-smad2/smad2  indicated  activation  of  related  signaling.  The 
immunobloting  has  three  replicates. 

ID2  functions  by  inhibiting  the  binding  to  DNA  of  E  proteins.  To  investigate  whether 
ID2  interacted  with  E  proteins  to  regulate  cell  fate  within  the  myeloid  lineage,  we 
analyzed  the  promoter  region  plus  the  5’  untranslated  region  of  KDR  for  putative  E 
protein  binding  sites  known  as  E  boxes.  We  identified  a  putative  binding  site  for  E2A 
(Figure  21a)  within  the  upstream  promoter  region  of  the  KDR  gene.  To  conclusively 
establish  the  ID2/E-protein/KDR  interaction,  we  performed  a  chromatin 
immunoprecipitation  assay  using  E2A  antibody  in  ID2-I-/-I-,  ID2  over-expressing,  and 
ID2-/-  bone  marrow  cells.  As  expected,  antibody  to  E2A  was  able  to  pull  down  the 
upstream  fragment  of  KDR.  Over-expression  of  ID2  abolished  the  binding  of  E2A  to  the 
promoter  of  the  KDR  gene  (Figure  21b).  Although  it  was  beyond  the  scope  and  focus  of 
our  analysis,  our  data  suggested  that  soluble  GM-CSF  and  TGF-P  may  influence  the 
expression  of  ID2  within  pluri-potential  bone  marrow  populations,  regulating  KDR 
expression  thereby  directing  cell  fate.  GM-CSF  and  TGF-P  induced  expression  of  ID2 
and  KDR  in  lineage  negative  bone  marrow  cells  (Figure  21c)  similar  to  conditioned 
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medium  derived  from  G1261  cells;  blocking  either  one  reversed  ID2  and  KDR 
expression.  These  results  support  the  hypothesis  that  glioma-secreted  GM-CSF  and  TGF- 
P  initiate  pro-tumoral  myeloid  cell  differentiation,  directing  the  formation  of  the  KDR-t- 
pro-angiogenic  immune  cell  population,  which  plays  a  crucial  role  in  the  tumor 
microenvironment  and  malignant  glioma  transformation. 
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4.  Key  Research  Accomplishments 

I.  We  have  demonstrated  that  myeloid  derived  suppressor  cells  increased  following  the 
progression  of  astrocytic  tumor  in  both  patients  and  murine  models. 

II.  We  identified  a  specific  population  across  endothelial  and  myeloid  lineages,  which  is 
VEGFR2+CD1  lb+  population  in  patients  and  tumor  bearing  mice. 

III.  We  successfully  performed  knock  out  VEGFR2  wihin  BMDCs  in  murine  glioma 
models,  and  observed  that  bone  marrow  derived  VEGFR2  contribute  to  tumor 
progression  and  animal  survival. 

IV.  We  have  performed  RNA-sequencing  on  tumor  associated  myeloid  progenitors,  and 
identified  inhibitor  of  DNA  binding  proteins  2  was  related  with  pro-tumoral  myeloid 
differentiation. 

V.  We  have  demonstrated  that  ID2  was  up-regulated  in  tumor  primed  hematopoietic 
progenitor  cells  in  both  patients’  samples  and  murine  models. 

VI.  We  identified  the  binding  region  of  E2A  in  VEGFR2  promoter,  which  indicates  role 
of  complex  E2A/ID2  in  regulation  of  VEGFR2  expression. 
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5.  Conclusion 


In  this  study,  we  have  utilized  glioma  patients  along  with  two  unique  murine  glioma 
models:  RCAS  glioma  model  and  G1261  model  to  study  the  BMDCs  during  different 
stages  of  glial  tumor.  Importantly,  we  identified  the  unique  the  population 
VEGFR2+MDSCs  in  both  patients  and  mice,  which  might  be  used  as  a  surrogate  marker 
for  glioma  diagnosis  and  prognosis  in  future.  We  have  validated  the  changes  of  myeloid 
lineage  and  endothelial  lineages  while  the  progression  of  gliomas,  and  observed  the 
increased  population  of  myeloid  derived  suppressor  cells  and  endothelial  progenitor  cells 
in  murine  glioma  models.  We  have  created  inducible  VEGFR2  knockout  system  in 
glioma  bearing  mice.  Taking  advantage  of  this  transgenetic  model,  we  demonstrated  that 
bone  marrow  derived  VEGFR2  signaling  plays  an  important  role  in  myeloid 
differentiation,  and  infiltration  into  tumor  tissues.  Deficiency  of  VEGFR2  in  BMDCs  led 
to  impairment  of  tumor  associated  myeloid  cells  and  delayed  progression  of  low-grade 
glioma.  All  of  these  findings  may  help  to  find  the  approach  to  suppress  the  progression  of 
low-grade  glioma  into  high-grade  form,  and  have  implications  to  predict  the  long-term 
survival  of  glioma  patients 

In  the  following  study,  we  will  keep  exploring  the  various  lineages  of  BMDCs  in  both 
low-grade  and  high-grade  glioma  patients.  We  would  further  validate  functional  role  of 
VEGFR2-I-  BMDCs  in  malignant  transformation,  and  investigate  the  connection  of 
myeloid  differentiation  with  tumor  associated  macrophages/neutrophils.  We  would 
dedicate  our  effort  to  delineate  the  signaling  pathways,  which  affect  pro-tumoral  myeloid 
cells,  and  further  characterize  the  downstream  signaling  of  inhibitor  of  DNA  binding 
proteins  2  m(ID2)  in  BMDCs  The  up-regulation  of  ID2  in  hematopoietic  progenitors 
provided  insight  how  tumor  enhances  VEGF  signaling  in  pro-tumoral  differentiation  via 
ID2/E2A  complex. 
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6.  Publications,  Abstracts,  and  Presentations 


1)  Abstract/Oral  presentation 

ID2/KDR  drives  the  differentiation  of  pro-malignant  myeloid  derived  suppressor  cells  in  glioma. 
Yujie  Huang,  Prajwal  Rajappa,  Jacqueline  Bromberg,  David  Lyden,  Jeffrey  Greenfield.  Cold 
Spring  Harbor-Asia  Meeting  (International)  “FRONTIERS  OF  IMMUNOLOGY  IN  HEALTH  & 
DISEASES”  September  2-September  6,  2014 

ID2/KDR  drives  formation  of  pro-malignant  myeloid  derived  suppressor  cells  in  glioma. 

Yujie  Huang,  Prajwal  Rajappa,  Jacqueline  Bromberg,  David  Lyden,  Jeffrey  Greenfield. 

Annual  Scientific  Meeting  of  the  Society  for  Neuro-Oncology,  November  13-16,  2014 

Oligodendrocyte  progenitor  cells  promote  neovascularization  in  glioma  by  disrupting  the  blood- 
brain  barrier.  Huang  Y,  Hoffman  C,  Rajappa  P,  Kim  JH,  Hu  W,  Huse  J,  Tang  Z,  Li  X,  Weksler  B, 
Bromberg  J,  Lyden  DC,  Greenfield  JP.  Cancer  Research.  2014  15;74(4):101 1-21. 


7.  Inventions,  Patents  and  Licenses 
None 


8.  Reportable  Outcomes 
None 


9.  Other  Achievements 
None 
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11.  Training  &  Professional  Development 


In  the  past  two  year,  I  have  received  extensive  trainings  and  related  proceedings  from 
various  workshops,  meetings,  and  hands-on  practices,  in  addition  to  regular  mentorships 
form  weekly  lab  meeting  and  journal  club.  To  follow  the  frontiers  of  tumor  immunology 
and  tumor  microenvironment  research,  I  have  attended  a  workshop  on  tumor 
microenvironment  (TME)  organized  by  national  cancer  institute  in  April  2014.  I  have 
communicated  several  renowned  experts  on  TME  with  their  study  and  our  proceedings.  I 
also  delivered  a  talk  about  our  work  in  Cold  Spring  Harbor-Asia  meeting  focusing  on 
immunology  in  diseases,  and  it  was  well  received.  In  November  2014,  my  abstract  on 
protumora  myeloid  differentiation  was  selected  for  oral  presentation  on  Annual  Meeting 
of  Society  of  Neuro-Oncology.  Part  of  my  research  on  PDGE  signaling  contributed  to  a 
publication  on  Cancer  Research  in  2014.  Additionally,  to  update  my  knowledge  on 
biomedical  and  genomic  fields,  I  have  continued  courses  and  workshops  including 
“Genomic  workshop”,  “Next  generation  sequencing  analysis”,  and  “Biostatics  for 
Clinical  Studies”,  which  are  provided  by  Clinical  Translational  Science  Center  in  Weill 
Cornell  Medical  College.  All  the  training  opportunities  armed  me  for  better  bench-side 
research  and  long-term  career  development.  According  to  the  feedback  from  my 
presentations  and  publications,  our  work  has  been  well  received  by  scientific  community. 
I  am  confident  to  achieve  the  career  goal  that  I  set  in  the  original  proposal. 
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Abstract 

Enhanced  platelet-derived  growth  factor  (PDGF)  signaling  in  glioma  drives  its  development  and  progression.  In 
this  study,  we  define  a  unique  role  for  stroma-derived  PDGF  signaling  in  maintaining  tumor  homeostasis  Avithin 
the  glioma  microenvironment.  Large  numbers  of  PDGF  receptor-a  (PDGFRa)-expressing  stromal  cells  derived 
from  oligodendrocytes  progenitor  cells  (OPC)  were  discovered  at  the  invasive  front  of  high-grade  gliomas,  in 
which  they  exhibited  a  unique  perivascular  distribution.  In  PDGFRa-deficient  host  mice,  in  which  orthotopic 
G1261  tumors  displayed  reduced  outgrowth,  we  found  that  tumor-associated  blood  vessels  displayed  smaller 
lumens  and  normalized  vascular  morphology,  with  tumors  in  host  animals  injected  with  the  vascular  imaging 
agent  gadolinium  also  being  enhanced  less  avidly  by  MRI.  Notably,  glioma-associated  OPC  promoted  endothelial 
sprouting  and  tubule  formation,  in  part  by  abrogating  the  inhibitory  effect  that  perivascular  astrocytes  exert  on 
vascular  endothelial  conjunctions.  Stromal-derived  PDGF-CC  was  crucial  for  the  recruitment  and  activation  of 
OPC,  insofar  as  mice  genetically  deficient  in  PDGF-CC  phenocopied  the  glioma/ vascular  defects  observed  in 
PDGFRa-deficient  mice.  Clinically,  we  showed  that  higher  levels  of  PDGF-CC  in  glioma  specimens  were 
associated  with  more  rapid  disease  recurrence  and  poorer  overall  survival.  Our  findings  define  a  PDGFRa/ 
PDGF-CC  signaling  axis  Avithin  the  glioma  stromal  microenAdronment  that  contributes  to  vascular  remodeling 
and  aberrant  tumor  angiogenesis  in  the  brain.  Cancer  Res;  74(4);  1011-21.  ©2013  AACR. 


Introduction 

The  tumor  microenAdronment  is  a  complex  amalgam  of 
tumor  cells,  intermingling  parenchymal  cells,  infiltrating 
immune  system  cells,  and  a  vascular  network  of  endothelial 
cells,  which  in  concert  determine  the  pathologic  and  biological 
features  of  distinct  tumors.  Cancer  diagnoses  and  grades 
significantly  vary  the  composition  of  the  tumor  microenAdron- 
ment  (1-3).  Gliomas,  malignances  arising  from  glial  cells  Avithin 
the  central  nervous  system  (CNS),  harbor  multiple  cell  types  in 
addition  to  cells  Avith  tumorigenicity.  During  normal  develop- 


Authors'  Affiliations:  ^Department  of  Neurological  Surgery,  The  Child¬ 
hood  Brain  Tumor  Project;  ^Department  of  Pediatrics;  ^Division  of  Hema¬ 
tology-Medical  Oncology;  "^Pediatric  Brain  Tumor  Research,  Children's 
Cancer  and  Blood  Foundation  Laboratories,  Weill  Cornell  Medical  College; 
Departments  of  ^Pathology  and  ^Medicine,  Memorial  Sloan-Kettering 
Cancer  Center,  New  York,  New  York;  and  ^National  Eye  Institute,  NIH, 
Bethesda,  Maryland 

Note:  Supplementary  data  for  this  article  are  available  at  Cancer  Research 
Online  (http://cancerres.aacrjournals.org/). 

Corresponding  Authors:  David  C.  Lyden,  51 5E  71st  Street,  New  York, 
NY  10021.  Phone:  212-746-3941;  Fax:  212-746-8423;  E-mail: 
dcl2001@med.cornell.edu;  and  Jeffrey  P.  Greenfield,  525E  68th  Street, 
New  York,  NY  10065.  Phone:  212-746-2363;  Fax:  212-746-7729;  E-mail: 
jpgreenf@med.cornell.edu 

doi:  10.11 58/0008-5472.CAN-1 3-1 072 

©2013  American  Association  for  Cancer  Research. 


ment  of  the  CNS,  proliferation,  differentiation,  and  migration  of 
glial  cells  are  tightly  controlled  and  well  described  (4,  5). 
However,  comparatively  little  is  knoAvn  about  the  composition, 
features,  and  functions  of  nontransformed  glial  cells  and  neural 
progenitor  cells  Avithin  the  natural  progression  of  glioma. 

The  involvement  of  platelet-derived  groAvth  factor  receptor- 
a  (PDGFRa)  in  gliomagenesis  has  been  well  demonstrated  (6). 
Upregulated  PDGF  signaling  through  PDGFRa  has  been  found 
to  be  a  common  feature  of  low-grade  astrocytic  and  oligoden- 
droglial  tumors  along  Avith  a  significant  subset  of  glioblastoma 
multiforme  (7,  8).  Interestingly,  the  adult  brain  contains  a 
Avidely  distributed,  abundant  progenitor  population  knoAvn  as 
oligodendrocyte  precursor  cells  (OPC),  which  have  been  sug¬ 
gested  to  have  high  tumor-initiating  potential  (9-12).  These 
cells  identified  by  expression  of  NG2/ 01ig2  are  normally  cycling 
and  express  PDGFRa  in  postnatal  brains  (13,  14).  Enhanced 
PDGF  signaling  is  present  in  a  specific  population  of  gliomas 
(15),  suggesting  that  OPCs  might  become  activated  and  under¬ 
go  specific  adaptations  in  response  to  unique  tumor  micro- 
environment  conditions.  In  this  study,  we  explore  the  distri¬ 
bution  and  possible  functional  contributions  of  OPCs  to  the 
vascular  microenvironment  dependent  on  PDGFRa/PDGF-CC 
signaling  during  glioma  progression.  We  suggest  a  role  for 
PDCFRa”^  OPCs  in  the  destabilization  of  vascular  architecture 
and  integrity  as  a  major  functional  contributor  to  tumor 
angiogenesis  and  tumor-grade  progression. 
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Materials  and  Methods 

Cells  lines  and  mice  lines 

deficient  mice  have  been  previously  described  (16). 
C57BL/6-derived  glioma  cells,  G1261  (authenticated  as  Supple¬ 
mentary  Table  SI)  and  KR158,  were  as  described  previously 
(17-19).  HCMEC/D3  have  been  previously  described  (20;  addi¬ 
tional  information  available  in  Supplementary  Methods). 

Intracranial  injections  and  tumor  imaging 

The  intracranial  injections  have  been  described  previously 
(21).  Bioluminesence  IVlS-100  (Xenogen)  and  MRl  (Bruker 
Biospin)  were  performed  to  monitor  the  progression  of 
tumor.  Tumor  margins  in  each  MRI  slice  were  manually 
outlined.  The  area  of  each  region  of  interest  was  calculated 
and  then  multiplied  by  the  slice  thickness.  All  slice  volumes 
were  added  up  to  calculate  the  volume  of  each  three-dimen¬ 
sional  (3D)  tumor.  For  detailed  methods,  please  refer  to 
Supplementary  Methods. 

Adoptive  bone  marrow  transplantation 

Transplantation  was  carried  out  as  previously  described  (22) 
Avith  slight  modifications.  Recipient  mice  were  lethally  irradi¬ 
ated  with  a  single  dose  of  9.5  Gy  to  their  whole  body  with  their 
heads  shielded.  Twenty-four  hours  after  irradiation,  5  x  10^ 
total  bone  marrow  cells  were  injected  by  tail  vein  followed  by  6 
weeks  to  allow  for  engraftment  of  bone  marrow  cells.  The 
engraftment  was  confirmed  by  complete  blood  count  and  fiow 
cytometry  to  ensure  adequate  hematopoiesis. 

Dextran  infusion  assay 

10  mg/mL  of  70-kD  rhodamine-conjugated  Dextran  (lysine 
fixable)  was  intravenously  injected  into  mice  (100  pL/mouse) 
as  previously  described  (23).  Two  hours  after  dextran  injection, 
mice  were  perfused  Avith  10  mL  PBS  followed  by  perfusion  of  5 
mL  of  4%  paraformaldehyde.  Brain  tissue  was  dissected  and 
fixed  for  the  subsequent  staining. 

Microarray  analysis 

Microarray  data  can  be  accessed  from  Gene  Expression 
Omnibus  (NCBI/GEO)  under  accession  number  GSE38283.  For 
detailed  procedures  of  the  preparation  of  samples,  please  refer 
to  Supplementary  Methods. 

Quantification 

All  data  are  given  as  mean  dz  SD  or  d=  SEM.  Differences  were 
compared  using  Student  ^  tests  or  one-way  ANOVA  followed  by 
post  hoc  tests.  For  more  experimental  procedures  and  quan¬ 
tification  details,  please  refer  to  Supplementary  Methods. 

Results 

stromal  cells  in  glioma  express  PDGFRa 

PDGF  signaling  is  upregulated  in  many  patients  Avith  glioma; 
30%  of  patients  harboring  high-grade  gliomas  express  an 
amplified  or  mutated  PDGFRa  gene  (7).  The  infiltration  of 
bone  marrow- derived  cells  (BMDC)  is  also  considered  as  a 
hallmark  of  tumor  progression  (22,  24).  To  better  understand 
the  temporal  and  spatial  variation  Avithin  the  glioma  micro- 


enAdronment  of  PDGFRa-expressing  cells  and  infiammatory 
BMDCs,  PDGFRa  staining  was  performed  on  orthotopic  G1261 
(murine-derived  glioma  cell  line)  tumors  in  mice,  which  had 
previously  been  adoptively  transplanted  Avith  GFP  bone  mar¬ 
row.  Large  numbers  of  BMDCs  and  PDGFRa^  cells  were 
identified  at  the  tumor  periphery  (Fig.  lA).  The  number  of 
these  peritumoral  PDGFRa^  cells  increased  at  day  14  along 
Avith  the  enhanced  infiltration  of  BMDCs  (Fig.  IB).  We  dem¬ 
onstrated  that  most  PDGFRa^  cells,  located  in  the  invasive 
front,  were  actually  nontumor  cells.  We  used  ortho  topic  glioma 
models  Avith  both  murine-derived  G1261  cells  (PDGFRa^) 
constitutively  expressing  GFP,  as  well  as  human-derived 
U87-MG  and  U251  cells  also  expressing  GFP  (Fig.  1C  and 
Supplementary  Fig.  SI  A),  to  show  that  these  were  not  invading 
tumor  cells.  Flow  cytometry  performed  upon  whole  tumor 
explants  confirmed  that  PDGFRa-expressing  cells  are  GFP 
negative  in  all  of  the  G1261,  U87-MG,  U251,  and  KR158  glioma 
models  (Fig.  ID  and  Supplementary  Fig.  SI).  The  existence  of 
peritumoral  PDGFRa-expressing  cells  was  also  verified  in  the 
KR158  glioma  model  (Supplementary  Fig.  SID).  Western  blot 
analyses  showed  enhanced  expression  of  PDGFRa  in  the 
tumor  margin,  compared  Avith  tissue  isolated  from  the  tumor 
core,  the  contralateral  hemisphere,  or  from  normal  age- 
matched  adult  mouse  brain  (Fig.  IE).  Immunostaining  of 
PDGFRa  showed  more  PDGFRa”^  cells  in  tumor  periphery 
than  in  the  contralateral  tissue  (Supplementary  Fig.  S2B). 
Taken  together,  these  results  demonstrate  that  large  numbers 
of  PDGFRa  stromal  cells  accumulate  in  the  periphery  and 
invade  the  front  of  gliomas. 

PDGFRa^  stromal  ceUs  derive  from  OPC  and 
predominate  perivascularly 

Different  tumors  are  supported  by  unique  microenviron¬ 
ments  containing  distinct  and  specific  stromal  cell  popula¬ 
tions.  To  better  define  the  origin  of  PDGFRa^  stromal  cells,  we 
studied  the  localization  of  PDGFRa  Avithin  G1261  tumor  speci¬ 
mens  using  other  lineage  markers  commonly  encountered 
Avithin  the  glioma  microenvironment.  PDGFRa^  stromal  cells 
did  not  express  markers  of  neural  stem  cells  (nestin),  astrocytes 
(glial  fibrillary  acidic  protein,  GFAP),  or  endothelium  (collagen 
IV),  nor  did  they  colocalize  Avith  BMDCs,  identified  using  GFP 
after  GFP^  bone  marrow  transplant  (Supplementary  Fig.  S2). 
PDGFRa^  cells  did,  however,  colocalize  Avith  NG2,  a  marker  for 
OPCs.  Moreover,  PDGFRa^  cells  display  ramified  processes, 
highly  suggestive  of  OPC  morphology,  and  expressed  oligoden- 
dritic  lineage  marker  01ig2  (Fig.  2A  and  Supplementary  Fig. 
S2C).  PDGFRa”^  glioma-associated  OPCs  (GA-OPC)  can  be 
identified  adjacent  to  vascular  endothelium  in  a  sandAvich-like 
configuration  (Fig.  2B).  Because  of  permeable  features  of  tumor 
endothelium  and  the  resolution  of  microscopy,  we  did  observe 
limited  overlaps  between  CD31 /PDGFRa.  The  GA-OPCs  do 
not,  however,  colocalize  with  pericytes  (PDGFRp^)  or  perivas¬ 
cular  astrocytes  (GFAP”^;  Fig.  2C  and  D).  In  addition  to  their 
abundance  Avithin  the  tumor  periphery,  OPCs  also  exist  Avithin 
tumor  cores  Avith  perivascular  processes  as  shoAvn  by  NG2- 
DsRed  (Supplementary  Fig.  S2D).  NG2  is  generally  used  as  OPC 
maker  in  CNS,  but  it  also  marks  certain  pericytes  and  microglial 
cells. 
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Figure  1 .  Stromal  cells  express  PDGFRa.  A,  PDGFRa  was  stained  in  GI261  tumors.  PDGFRa^  cells  (red)  and  GFP*^  BMDCs  (green)  are  inset  (left,  tumor  center; 
right,  tumor  periphery).  Scale  bar,  1 00  |j,m,  n  =  4  mice.  B,  GI261  tumors  in  C57BL/6  mice  (GFP^  BMDCs)  were  harvested  at  D1 ,  D7,  D1 4,  and  D21  postinjection. 
The  number  of  PDGFRa  cells  and  GFP^  BMDCs  in  the  tumor  periphery  was  quantified  (cell  count  per  field).  Means  ±  SD,  n  =  5  mice.  Cell  counts  in  D1 4  and 
D21  were  significantly  higher  than  those  in  D1  and  D7.P<  0.001  by  one-way  ANOVA.  C,  PDGFRa  was  stained  (red)  in  GFP-GL261  (green)  tumors  (day  21). 
Tumor  periphery  (left)  and  invasive  tumor  margin  (right)  are  shown.  Scale  bar,  50  |j,m.  D,  flow  cytometry  was  performed  on  dissociated  cells  from 
GFP-GI261  tumors  with  isotype  IgG  (left)  and  PDGFRa  antibody  (right).  E,  tissues  from  tumor  core,  tumor  margin,  contralateral  hemisphere,  and  normal  adult 
murine  brain  were  blotted  by  PDGFRa  and  phos-PDGFRa.  GFP  and  GAPDH  were  used  as  loading  controls  with  tumor  cells  expressing  GFP. 

Normalized  densitometric  analyses  on  three  independent  experiments  (n  =  3  mice)  were  provided  below  each  lane.  The  level  of  PDGFRa  in  tumor  margin 
is  significantly  higher  than  those  in  other  groups  (P  <  0.001)  by  one-way  ANOVA. 


PDGFRa^  stromal  cells  facilitate  the  progression  and 
neovascularization  of  glioma 

Tumor  cells  may  educate  adjacent  nontransformed  cells 
into  forming  tumor-associated  stroma,  which  can  be  differen¬ 
tiated  from  the  normal  tissue  (3,  25).  In  response  to  malig¬ 
nancy,  OPCs  may  be  activated  through  PDGFRa  signaling 
(9,  26).  To  understand  the  functional  role  of  GA-OPCs  in  the 
progression  of  glial  malignancies,  we  created  an  orthotopic  and 
syngeneic  model  using  G1261  cells  transfected  with  a  luciferase 
gene.  These  cells  were  implanted  in  PDGFRa-knockout  mice  to 
study  the  importance  of  PDGFRa  signaling  between  stromal 
and  tumors  cell  populations,  upon  glioma  progression.  In  Rosa- 
cre/PDGFRo/^^^  mice,  growth  of  G1261  tumors  was  inhibited 
compared  with  tumor  growth  in  Rosa-cre/PDGFRa^^^  mice 
(Fig.  3A).  MRI  analysis  of  G1261  tumors  in  PDGFRc/^^  mice 
demonstrated  animals  with  less  tumor  burden  and  tumors 


with  less  contrast  enhancement  on  MRI  than  tumors  within 
control  mice  (Fig.  3B).  Compared  with  tumors  in  PDGFRa^^^ 
mice,  tumors  in  PDGFRo/^^^  mice  had  a  slightly  higher  blood 
vessel  density  (Supplementary  Fig.  S3A)  and  demonstrated 
distinct  vascular  morphology.  Gliomas  within  PDGFRc/^^^  mice 
demonstrated  considerably  smaller  diameters  of  vessel  lumens 
(Fig.  3C),  similar  to  vasculature  typical  within  the  normal 
mouse  brain.  Importantly,  within  PDGFRc/^^^  mouse  tumors, 
pericyte  density  and  architecture  was  not  significantly  changed 
from  PDGFRa^^^  mice  as  visualized  by  the  staining  for 
PDGFRp  (Supplementary  Fig.  S3B). 

Neovascularization  of  gliomas  is  distinct  from  other  types  of 
tumors.  Perhaps  due  to  the  brain's  robust  vascular  network, 
brain  tumors  either  remodel  preexisting  blood  vessels  or 
recruit  endothelial  progenitor  cells  to  undergo  microvascular 
proliferation,  a  hallmark  of  high-grade  gliomas  (1,  27). 
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Figure  2.  OPC-derived  PDGFRa 
stromal  cells  are  perivascular.  A, 
PDGFRoc  (red)  colocalizes  with 
NG2  (green)  in  the  tumor  periphery. 
Arrow,  overlap.  Asterisk,  NG2^ 
cells  (green)  with  suggestive 
pericyte  morphology.  Scale  bar, 
50  |j,m.  B,  PDGFRa  is  adjacent  to 
blood  vessel  marker  CD31  and  is 
not  colocalized  with  the  pericyte 
marker  PDGFRp  (C).  Z-stack 
imaging  was  used  with  confocal 
microscopy.  Scale  bar,  20  |a,m.  D, 
costaining  of  GFAP  (perivascular 
astrocytes)  and  PDGFR(3 
(pericytes)  on  GI261  tumor 
periphery.  Arrows,  vascular 
lumens;  scale  bar,  20  |im. 


Compared  with  normal  brain  and  low-grade  gliomas,  grade  III 
and  IV  gliomas  have  distorted,  enlarged  blood  vessels  but  do 
not  necessarily  demonstrate  higher  blood  vessel  density  (27). 
Blood  vessels  in  high-grade  gliomas  are  distinguishable  from 
those  in  low-grade  gliomas  not  by  their  density,  but  by  their 
morphology.  In  PDGFRc/’^^^  mice,  G1261  brain  tumor-associ¬ 
ated  blood  vessels  were  not  enlarged  or  distorted  compared 
with  that  within  tumors  in  wild-type  mice.  Tumor  cells  in 
PDGFRc/^^^  mice  have  less  mitotic  activity  (decreased  Ki67 
staining),  less  hypoxia  demonstrated  (decreased  pimonidazole 
staining),  and  less  leakiness  as  indicated  by  a  dextran  infusion 
assay  (Supplementary  Fig.  S3C).  Other  factors,  which  could 
theoretically  induce  changes  in  vascular  morphology,  were 
studied  in  both  groups  of  animals  by  measuring  the  infiltration 
of  BMDCs  and  monitoring  hematopoiesis;  there  were  moder¬ 


ate  decreases  of  infiltrating  BMDCs  within  G1261  tumors  in 
PDGFRc/^^^  mice  (Supplementary  Fig.  S4A),  and  hematopoiesis 
does  not  appear  to  be  affected  (Supplementary  Fig.  S4B). 

GA-OPCs  facilitate  angiogenesis 

To  better  understand  the  contribution  of  GA-OPCs  to  the 
vascular  remodeling  occurring  during  glioma  neovasculariza¬ 
tion,  we  performed  in  vitro  assays  using  primary  GA-OPCs 
cultures.  Primary  human-derived  GA-OPCs  and  glioma  cells 
were  isolated  directly  from  glioma  resections  (Supplementary 
Fig.  S5A);  OPCs  were  maintained  in  the  progenitor  stage  within 
the  oligodendritic  lineage  (Supplementary  Fig.  S5B  and  S5C). 
Using  a  3D  culture  system,  we  cocultured  human  brain  endo¬ 
thelial  cells  (HCMEC/D3)  with  GA-OPCs  or  human  astrocytes, 
or  both,  to  study  their  interactions.  We  labeled  the  HCMEC/D3 
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Figure  3.  Knocking  out  PDGFRa  in  stromai  ceiis  inhibits  giioma  progression.  Luciferase-Gi261  ceiis  were  intracraniaiiy  injected  into  RosaCre-ERT2/ 
PDGFRa^'^^'  mice  and  RosaCre-ERT2/PDGFRa‘^^‘^  mice.  A,  tumor  burden  is  represented  by  photon  fiux  intensity.  Means  ±  SD,  n  =  1 0  mice;  f  test,  *,P<  0.05. 
B,  precontrast  and  postcontrast  T1 -weighted  MRis  were  performed  on  each  group  at  day  15.  Arrows,  enhancing  tumor.  Tumor  voiume  quantification  is 
demonstrated.  Means  ±  SD,  n  =  5  mice;  f  test,  **,  P<  0.01 .  C,  vWF  was  stained  in  the  tumor  periphery  and  within  the  tumor  of  each  group  (ieft  and  middie 
paneis).  Coiiagen  iV  was  stained  in  each  group  (right).  Scaie  bar,  50  pm.  Vessei  diameters  were  quantified.  Means  ±  SD,  n  =  6  mice;  t  test,  *,  P  <  0.05. 


cells  with  GFP  and  allowed  them  to  form  tubules.  Once  a  tubule 
network  was  established,  equal  numbers  of  GA-OPCs  and 
astrocytes  were  introduced  to  the  endothelial  cell  cultures. 
The  in  vitro  3D  coculture  system  recapitulated  the  in  vivo 
phenomenon  within  which  the  astrocytic  foot  processes  sur¬ 
round  endothelium  to  form  the  blood-brain  barrier.  We  found 
that  both  astrocytes  and  GA-OPCs  were  adjacent  to  the 
endothelial  tubules.  Astrocytic  foot  processes  formed  an  intact 
and  continuous  tubule-like  structure  close  to  endothelial 
tubules.  In  contrast,  GA-OPCs  appeared  to  integrate  with  the 
endothelial  tubules  but  did  not  form  continuous  intercellular 
interaction  with  each  other  as  did  astrocytes  (Fig.  4A).  Inter¬ 
estingly,  if  we  allowed  astrocytes  to  establish  interaction  with 
HCMEC/D3  first  with  coculture,  and  then  added  GA-OPCs,  the 
GA-OPCs  could  disrupt  the  interaction  formed  between  astro¬ 
cytes  (Fig.  4A).  We  further  used  in  vitro  models  to  test  the 
permeability  of  these  intercellular  interactions,  which  were 
formed  by  endothelial  cells  under  differing  coculture  condi¬ 
tions.  These  assays  demonstrated  that  astrocytes  protect  the 
integrity  of  the  endothelial  cell  monolayer  while  bathed  in  U251 
cell-derived  conditioned  medium,  whereas  GA-OPCs  abrogate 
the  protective  effect  of  astrocytes  by  creating  a  more  perme¬ 
able  endothelial  cell  monolayer  (Supplementary  Fig.  S6C). 

In  addition  to  defining  interactions  between  endothelial 
cells,  astrocytes,  and  GA-OPCs,  we  performed  a  series  of 
angiogenesis  assays  to  study  the  mechanism  through  which 
GA-OPCs  impact  vascular  remodeling.  We  used  a  sprouting 
assay  to  investigate  the  infiuence  of  perivascular  stromal  cells 
upon  the  migration  of  endothelial  cells.  Astrocytes  suppressed 


sprouting  from  endothelial  spheres  made  from  HCMEC/D3; 
GA-OPCs  reversed  this  inhibitory  effect  when  cocultured  with 
astrocytes  (Fig.  4B).  We  also  analyzed  the  number  of  endothe¬ 
lial  cells  in  coculture  with  astrocytes  and  GA-OPCs.  Both 
astrocytes  and  GA-OPCs  only  minimally  affected  the  number 
of  endothelial  cells  (Supplementary  Fig.  S5D),  suggesting  that 
the  regulatory  effect  of  stromal  cells  on  endothelial  sprouting 
was  proliferation  independent.  This  function  is  PDGFRa 
dependent;  a  neutralizing  PDGFRa-specific  antibody  reversed 
the  effect  (Supplementary  Fig.  S6A).  Tubule  formation  assays 
similarly  demonstrated  that  astrocytes  inhibited  the  formation 
of  new  endothelial  tubules  and  GA-OPCs  reversed  this  in  a 
PDGFRa  signaling-dependent  manner  (Fig.  4C  and  Supple¬ 
mentary  Fig.  S6B). 

How  these  interactions  between  endothelial  cells  and  OPCs 
directly  affect  the  tumor  progression  is  difficult  to  precisely 
define.  We  have  demonstrated  that  inhibiting  OPCs  clearly 
affects  angiogenesis,  which  might  then  inhibit  tumor  progres¬ 
sion.  More  interestingly,  there  may  also  exist  a  direct  link 
between  endothelial  and  tumor  cells.  We  hypothesized  that 
OPCs  may  direct  endothelial  cells  into  a  more  active  state 
through  upregulation  of  angiocrine  factors,  thereby  promoting 
malignancy.  We  cocultured  endothelial  cells  with  G1261  cells  in 
the  presence  or  absence  of  VEGF  (to  stimulate  endothelial 
cells).  We  found  that  tumor  cells  proliferated  more  when 
cocultured  Avith  activated  endothelial  cells  than  when  cocul¬ 
tured  with  normal  endothelial  cells  (Supplementary  Fig.  S7). 
Considering  the  complex  nature  of  the  glioma  microenviron¬ 
ment,  interaction  between  OPCs  and  endothelial  cells  may  be 
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Figure  4.  GA-OPCs  promote  opening  of  the  endothelial-astrocytic  junction  and  angiogenic  activity  of  endothelial  cells.  A,  GFP^  HCMEC/D3  cells  (green)  were 
seeded  in  Matrigel  to  assess  tubule  formation.  Rodamine  labeled  (red),  human  astrocytes  (HA),  and  GA-OPC  were  added  subsequently.  For  the  combination  of 
human  astrocytes  and  OPCs,  labeled  human  astrocytes  were  allowed  to  seed  with  endothelial  tubules  first,  followed  by  nonlabeled  GA-OPCs  12  hours  later. 
Scale  bar  of  the  first  left  panel,  1 00  Scale  bar  of  right  three  panels,  50  |j,m.  GFP^  endothelial  tubules'  length  and  coverage  were  quantified.  Means  ±  SEM, 
n  =  5;  one-way  ANOVA,  *,  P  <  0.05.  B,  coculture  with  human  astrocytes,  GA-OPC  (glioma  associated  OPCs),  or  both.  The  coculture  system  was 
maintained  in  U251 -conditioned  medium.  Endothelial  sprouting  was  quantified  by  the  area  of  endothelial  cells  beyond  the  primary  spheroid  (dashed  lines). 
Means  ±  SEM,  n  =  6;  one-way  ANOVA,  ***,  P<  0.001 . 0,  GFP-HCMEC/D3  cells  were  seeded  on  Matrigel  to  form  tubules,  cocultured  with  GA-OPC,  human 
astrocytes,  or  both,  and  maintained  in  U251  -conditioned  medium.  The  GFP"^  endothelial  tubules  in  each  group  were  quantified  by  length.  Means  ±  SEM,  n  =  6; 
one-way  ANOVA,  **,P<  0.01 .  All  graphs  were  acquired  from  two  independent  experiments.  N,  the  number  of  independent  wells  of  treatment  per  condition. 


even  more  complicated  than  that  hypothesized  with  multiple 
pathways  independently  promoting  malignancy. 

BMDCs  and  stromal-derived  PDGF-CC  are  key  mediators 
to  GA-OPCs 

In  the  context  of  the  glioma  microenvironment,  OPCs 
exponentially  increase  in  number  and  in  the  expression  of 
PDGFRa  (Fig.  1).  Signaling  that  originates  from  tumor- 
initiating  or  alternative  stromal  cells  might  trigger  OPC 


activation  and  redirect  them  into  protumoral  phenotypes. 
As  described  above,  the  number  of  PDGFRa^  cells  dramat¬ 
ically  increased  around  day  14  after  injection  of  tumor  cells, 
which  mirrored  the  pattern  of  immune  cell  infiltration  (Fig. 
lA  and  B).  Therefore,  we  sought  to  test  the  hypothesis  that 
infiltrating  BMDCs  contribute  to  the  recruitment  and  acti¬ 
vation  of  GA-OPCs  during  tumor  progression.  In  addition, 
BMDCs  are  detected  not  only  intratumorally,  but  also  within 
the  tumor  periphery  (Fig.  lA),  a  localization  that  would 
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Figure  5.  Stromal  derived  PDGF-CC  serves  as  an  activator  of  OPCs  in  glioma.  A,  PDGF-C  is  upregulated  in  glioma-associated  macrophages.  Quantitative  RT- 
PCR  of  PDGF-C  on  GFP+/F4/80+  macrophage/microglia  sorted  from  orthotopic  GI261  tumors  or  normal  brain.  Means  ±  SEM,  n  =  6  mice;  f  test,  ***,  P  < 
0.001 .  B,  the  growth  of  GI261  tumors  in  mice  transplanted  with  POGF-C”"^”  bone  marrow  or  littermate  PDGF-C"'''^''’  bone  marrow  was  measured  by  photon  flux. 
Means  ±  SD,  n  =  8  mice;  t  test,  * ,  P  <  0.05.  C,  PDGFRoc-f  cells  were  quantified  in  the  tumor  periphery  of  GI261  tumors  at  D1 4  and  D21 ,  respectively,  from  both 
PDGF-C“^“  and  PDGF-C^^^  BMT  animals.  Means  ±  SEM,  n  =  5  mice;  one-way  ANOVA,  **,  P  <  0.01.  D,  endothelial  cells  and  astrocytes  isolated  from 
cortical  tissue  of  C57BL/6  mice  were  treated  with  GI261  conditioned  medium  or  control  medium.  Quantitative  RT-PCR  of  PDGF-C  was  performed, 
demonstrating  the  upregulation  of  PDGF-C  in  both  endothelial  and  astrocytic  cells  in  response  to  glioma  conditions.  Means  ±  SEM,  n  =  6;  f  test,  *** ,  P<  0.001 . 
The  graph  was  acquired  from  two  independent  experiments.  N,  the  number  of  independent  wells  of  treatment  per  condition. 


permit  direct  BMDCs/GA-OPCs  interactions.  Subsequently, 
we  further  characterized  the  inhitrating  GFP^  BMDCs  in 
G1261 -implanted  mice,  and  found  that  more  than  80%  of  the 
BMDCs  in  these  tumors  are  Fd/SO”^  macrophages/microglial 
cells  (Supplementary  Fig.  S8A).  We  then  performed  micro¬ 
array  analysis  on  GFP'^/F4/80'^  BMDCs,  and  found  that  the 
PDGF-C  gene  was  signihcantly  upregulated  in  glioma-asso¬ 
ciated  BMDCs  compared  with  BMDCs  from  normal  murine 
brain  (Supplementary  Fig.  S8B).  PDGF-C  is  a  cytokine  that 
preferentially  homodimerizes  into  PDGF-CC,  which  is  the 
ligand  predominantly  binding  PDGFRa.  The  upregulation  of 
PDGF-C  in  bone  marrow-derived  macrophage /microglia 
isolated  from  murine  gliomas  was  conhrmed  by  real-time 
PCR  (RT-PCR;  Fig.  5A). 

PDGF-C  deletion  recapitulates  PDGFRa  null  phenotype 

To  examine  whether  PDGF-CC  derived  from  bone  marrow- 
derived  macrophages  or  microglial  cells  might  push  GA-OPCs 
into  tumor-supporting  phenotyes,  we  transplanted  wild  type 
C57BL/ 6  mice  with  bone  marrow  from  PDGF-C~^~  donor  mice. 
Six  weeks  after  bone  marrow  engraftment,  luciferase-labeled 
G1261  cells  were  injected  intracranially  into  recipient  mice. 
PDGF-C~^~  mice  did  not  show  any  developmental  defects  in 


their  hematopoietic  system  (Supplementary  Fig.  S9);  however, 
we  found  that  tumors  implanted  within  mice  transplanted 
with  PDGF-C~^~  bone  marrow  grew  more  slowly  at  an  early 
stage  of  tumorigenesis  (day  14)  than  tumors  in  mice  trans¬ 
planted  with  littermate  PDGF-C^^^  bone  marrow.  However,  in 
later  stages  (day  21)  of  tumor  growth,  PDGF-C~^~  bone  mar¬ 
row-transplanted  animals  no  longer  demonstrated  a  statisti¬ 
cally  signihcant  reduction  of  tumor  size  (Fig.  5B).  Within  the 
tumor  periphery  specihcally,  we  found  signihcantly  fewer 
PDGFRa^  cells  in  the  PDGF-Cr^~  bone  marrow-transplanted 
group  at  D14  (P  =  0.0051),  but  not  at  a  later  stage  (D21;  P  = 
0.32;  Fig.  5C).  These  results  suggest  that  microphage-/micro- 
glia-derived  PDGF-CC  did  contribute  to  the  recruitment  and 
activation  of  GA-OPCs,  a  process  not  necessarily  dependent 
upon  bone  marrow-derived  macrophages/microglia  at  later 
stages  of  tumor  development. 

Soluble  PDGF-CC  may  not  be  solely  secreted  by  BMDCs; 
thus,  we  inferred  that  other  stromal  cells  within  the  glioma 
microenvironment  such  as  residential  microglia,  activated 
astrocytes,  or  angiogenic  endothelial  cells  could  be  candidates 
as  sources  of  soluble  PDGF-CC.  Therefore,  we  tested  PDGF-C 
expression  in  astrocytes  and  endothelial  cells  in  response  to 
tumor  cells.  Primary  cultures  of  astrocytes  and  endothelial 
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Figure  6.  Glioma  growth  in  PDGF-C  deficient  mice  is  phenotypically  similar  to  glioma  growth  in  PDGFRa  KO  mice.  A,  growth  of  orthotopic  GI261  tumors  in 
PDGF-C“^“  and  PDGF-C^^^  mice.  T umor  burden  is  represented  by  bioluminescence  photon  flux.  Means  ±  SD,  n  =  8  mice;  f  test,  ** ,  P<  0.01 .  B,  pre-  and  post¬ 
contrast  T1 -weighted  MRI  was  used  to  image  GI261  tumors  in  PDGF-C“'^“  animals  and  PDGF-C'^'''^  animals  on  day  17.  Arrows,  the  enhancing 
tumor  region.  Tumor  volumes  were  quantified.  Means  ±  SD,  n  =  5  mice;  f  test,  **,  P  <  0.01 .  C,  PDGFRa”^  cells  were  decreased  in  PDGF-C“'^“  mice.  Blood 
vessels  (collagen  IV)  and  PDGFRa  were  stained  in  both  the  invasive  tumor  front  (invading  front)  and  from  within  GI261  tumors  (intratumor)  grown  in  the 
PDGF-C“^“  group  and  the  PDGF-C^^^  group.  Scale  bar,  50  |im.  Means  ±  SD,  n  =  5  mice;  t  test,  p  <  o.OOl . 


cells  were  isolated  from  mouse  brain,  and  treated  with  or 
without  conditioned  medium  derived  from  G1261  cell  culture. 
RT-PCR  from  both  primary  cultures  showed  that  the  expres¬ 
sion  of  PDGF-C  was  upregulated  by  tumor-conditioned  medi¬ 
um  (Fig.  5D),  implying  that  PDGF-CC  might  be  expressed  and 
secreted  by  glioma-associated  astrocytes  and  endothelial  cells. 
However,  we  found  that  other  types  of  tumor  cells,  including 
B16/F10  and  KR158,  could  both  upregulate  the  PDGF-C  expres¬ 
sion  within  stromal  cells  (Supplementary  Fig.  SlOA),  suggest¬ 
ing  that  the  upregulation  of  PDGF-C  may  be  specific  to  glioma 
stroma,  but  not  to  tumor  cells.  The  stromal-derived  PDGF-CC 
may  also  be  related  to  other  cancers  and  proinfiammatory 
states. 

To  assess  the  contribution  of  total  stromal  PGDF-CC  to 
gliomagenesis,  we  intracranially  injected  G1261  cells  into 
PDGF-C~^~  knockout  mice.  G1261  tumor  growth  in  PDGF- 
C~^~  mice  was  significantly  slower,  and  tumor  burden  was 
less  than  the  tumors  in  littermate  PDGF-C^ mice  (Fig.  6A  and 
B).  MRI  demonstrated  strikingly  less  tumor  enhancement  in 
PDGF-C~^~  mice.  These  results  recapitulated  the  phenotype  of 
G1261  tumor  progression  in  PDGFRa  knockout  mice.  More 
importantly,  we  found  many  fewer  PDGFRa”^  cells  in  the  tumor 
periphery  in  PDGF-Cr^~  mice  and  significantly  smaller  vascu¬ 
lar  lumens  (Fig.  6C).  Taken  together,  these  results  demonstrate 


that  stromal-derived  PDGF-CC  is  a  key  mediator  in  the  recruit¬ 
ment  and  activation  of  OPCs  during  the  progression  of  gliomas. 

PDGF-CC,  of  which  the  primary  receptor  is  PDGFRa,  has 
been  implicated  in  various  cell  types  and  pathologic  condi¬ 
tions.  A  recent  study  showed  that  PDGF-CC  knockout  mice 
showed  abnormal  cerebral  vascularization  (28),  supporting  the 
role  of  PDGF-CC  in  regulating  the  vascular  architecture  under 
physiologic  conditions  in  the  CNS.  In  the  postnatal  brain, 
PDGFRa  is  exclusively  expressed  on  OPCs  (29).  Thus,  the 
stromal  PDGF-CC/PDGFRa  axis  maybe  an  important  intrinsic 
angiogenic  signaling  pathway  regulating  angiogenesis  in  glio¬ 
mas— independent  of  the  classic  VEGF-mediated  pathway.  To 
further  understand  the  role  of  this  signaling  axis  in  regulating 
angiogenesis,  we  used  an  angiogenesis  antibody  array  (R&D 
systems)  to  profile  the  angiogenic  factors  derived  from  OPCs 
after  the  exposure  to  PDGF-CC.  PDGF-CC  upregulated  proan- 
giogenic  factors  such  as  CXCLIO,  CX3CL1,  and  IGFBP-1  in 
OPCs  (Supplementary  Fig.  SlOB).  PDGF-AA,  PDGF-AB,  PDGF- 
BB,  and  VEGF  remained  unchanged  in  the  same  assay.  Among 
upregulated  factors,  OPC-derived  MMP9  was  recently  found  to 
mediate  blood-brain  barrier  opening  in  response  to  white 
matter  injury  (26).  These  results  indicate  that  GA-OPCs  might 
have  a  direct  impact  upon  activating  endothelial  cells  and 
recruiting  proangiogenic  myeloid  cells. 
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Figure  7.  Human  glioblastoma 
(GBM)  samples  from  TCGA 
support  a  role  for  the  PDGF- 
signaling  axis.  A,  patient  samples 
show  GA-OPCs.  Sections  from 
IDH  mutant-positive/grade  III 
astrocytomas  were  stained  with  an 
antibody  against  IDH  R132H  and 
PDGFRa.  Arrow  demonstrates 
PDGFRa+  IDH  R132H  cells.  The 
inset  shows  an  enlarged  image. 
Scale  bar,  50  |j,m.  B,  PDGF-C  levels 
are  reversely  correlated  with 
glioblastoma  patient's  survival. 
The  overall  survival  and 
progression-free  survival  of 
patients  with  glioblastoma  were 
analyzed  based  on  the  expression 
of  PDGF-C  by  using  data  from 
TCGA  in  the  classic  subtype  of 
glioblastoma.  Survival,  P  =  0.008; 
progression-free  survival,  P  = 
0.006  by  the  log-rank  test. 
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GA-OPCs  in  patients  with  glioma 

The  brain  tumor  microenvironment  in  patients  with  glioma 
is  very  complex.  The  signaling  within  undefined  stromal 
populations  and  tumor  cells  is  intricate.  This  is  further 
complicated  by  the  molecular  heterogeneity  between  subpo¬ 
pulations  of  actual  glial  tumor  cells.  Therefore,  we  sought  to 
confirm  and  quantify  the  existence  of  PDGFRa”^  stromal  cells 
directly  from  samples  resected  from  patients  diagnosed  with 
a  high-grade  glioma.  Among  established  markers  for  glioma 
cells,  isocitrate  dehydrogenase  1  (IDHl)  mutation  is  extreme¬ 
ly  reliable  in  distinguishing  tumor  cells  from  stromal  cells  as  it 
has  been  used  to  detect  even  single  disseminated  tumor  cells 
(30,  31).  Grade  III  astrocytomas  bearing  IDHl  mutations  were 
chosen,  and  stromal  cells  were  characterized  as  negative  for 
the  IDH  mutation.  Doubly  staining  for  PDGFRa  and  IDH 
(R231H),  we  demonstrated  a  population  of  cells  positive  for 
PDGFRa  and  negative  for  IDH  (R231H;  Fig.  7A).  Interestingly, 
this  region  was  actively  undergoing  microvascular  prolifera¬ 
tion  evidenced  by  enlarged  vascular  lumens  (Supplementary 
Fig.  Sll). 

As  we  demonstrated  in  our  animal  glioma  models,  PDGF-CC 
is  a  potential  ligand  through  which  the  activation  of  GA-OPCs 
may  occur,  thereby  contributing  to  angiogenesis  and  glioma 
progression.  In  PDGF-C  null  mice,  glioma  progression  was 
slowed  due  to  a  lack  of  activated  GA-OPCs.  We  examined  the 
expression  of  PDGF-C  in  different  types  of  gliomas  within 
Rembrandt.  We  observed  that  expression  levels  of  PDGF-C 
were  higher  in  high-grade  gliomas  as  compared  with  low-grade 
gliomas  or  nontumor  diseases  (Supplementary  Fig.  SI2).  To 
further  study  the  role  of  PDGF-CC  in  human  high-grade  glio¬ 
mas,  data  from  The  Cancer  Genome  Atlas  (TCGA)  for  different 


subsets  of  patients  with  glioblastoma  were  analyzed  based  on 
their  expression  of  PDGF-C.  Samples  were  subgrouped  into 
those  with  higher  PDGF-C  expression  (top  10-14%)  with  the 
remainder  classified  as  PDGF-C  low.  Analyses  of  survival  and 
progression  were  performed  within  four  recognized  subsets  of 
glioblastoma,  including  classic,  proneural,  mesenchymal,  and 
neural.  Within  the  classic  subset  of  glioblastoma,  but  not  in 
other  subsets,  expression  of  PDGF-C  was  significantly  correlat¬ 
ed  with  disease  progression  {P  =  0.0060),  and  inversely  corre¬ 
lated  with  survival  of  patients  {P  =  0.0081;  Fig.  7B  and  Sup¬ 
plementary  Fig.  SI3).  This  demonstration  within  the  classic 
subset  of  human  glioblastoma  is  consistent  with  our  demon¬ 
strated  slowing  of  orthotopic  G126I  tumor  progression  in 
PDGD-C  null  mice.  Surprisingly,  even  other  members  of  the 
PDGF  family,  including  PDGF-A  and  PDGF-B,  commonly  upre- 
gulated  in  glioblastoma,  were  not  found  to  be  statistically 
significantly  {P  >  0.05)  correlated  with  patients'  survival  or 
disease  progression  (Supplementary  Fig.  SI4). 

Discussion 

OPCs  are  the  most  abundant  neural  progenitor  cells  in 
postnatal  brains.  They  can  differentiate  into  oligodendrocytes 
or  astrocytes  in  response  to  different  stimuli  to  aid  in  myelina- 
tion  and  wound  healing,  among  other  roles  in  the  CNS  (32-36). 
Recent  studies  also  suggest  that  OPCs  have  a  higher  degree  of 
plasticity  and  may  be  more  sensitive  to  transformation  (9, 10). 
In  certain  rodent  brain  tumor  models,  it  has  been  demonstrat¬ 
ed  that  recruited  stromal  OPCs  can  be  tumorigenic.  However, 
their  transformation  of  OPCs  may  require  a  preexisting  onco¬ 
genic  mutation  such  as  Ink4a/Arf  or  P53  as  indicated  by  these 
studies  (37). 
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Invasive  glioma  cells  frequently  migrate  along  myelinated 
white  matter  fiber  tracts  (38),  a  phenomenon  whose  underlying 
mechanism  is  not  clearly  elucidated.  It  has  also  been  estab¬ 
lished  that  OPCs  preferentially  reside  along  identical  pathways 
of  dissemination  (39,  40).  Many  studies  have  examined  and 
demonstrated  that  glioma  progression  and  invasion  follows  a 
perivascular  pattern  (25,  41,  42),  particularly  at  disease  recur¬ 
rence.  In  this  study,  we  demonstrated  angiogenic  endothelial 
cells,  perivascular  macrophages /microglia,  and  OPCs  working 
in  concert  to  form  a  proangiogenic  and  proinvasive  niche  at  the 
invading  front  of  glioma.  Our  data  suggest  that  stromal  OPCs 
are  a  key  element  not  only  in  initiating  angiogenesis  but  also  in 
driving  glioma  invasion.  To  extend  our  hypotheses  derived 
from  mouse  models  to  humans,  we  selected  IDH1R132H 
mutant  grade  III  astrocytomas  Avithin  which  to  analyze  stromal 
PDGFRa^  cells  from  human  specimens.  The  IDHl  R132H 
mutation  is  an  ideal  marker  to  distinguish  between  tumor 
cells  and  stromal  cells  in  grade  III  astrocytomas.  In  addition, 
grade  III  astrocytomas  are  by  definition  in  the  transition  stage, 
during  which  we  suggest  that  OPCs  play  a  crucial  role  in  both 
microvascular  proliferation  and  tumor  invasion. 

In  the  previous  work,  PDGF-CC  was  found  to  induce  the  loss 
of  the  blood-brain  barrier  and  induce  vascular  leakiness  in  a 
stroke  model,  and  PDGF-C  deficiency  leads  to  abnormal 
cerebral  vascularization  (28, 43).  In  our  study  in  tumor-bearing 
mice,  the  number  of  OPCs  Avith  activated  PDGFRa  was 
increased  compared  Avith  non-tumor-bearing  mice.  However, 
similar  tumors  could  not  increase  the  number  of  OPCs  in 
PDGF-C“^“  mice.  Interestingly,  PDGF-CC  has  been  identified 
as  a  key  factor  alloAving  tumors  to  recur  folloAving  the  anti- 
VEGF  treatment  in  glioblastoma  and  other  types  of  cancer  (44, 
45).  All  of  this  eAddence  suggests  that  PDGF-CC  plays  an 
important  role  in  the  vascular  remodeling  occurring  in  glio¬ 
mas,  which  may  be  independent  of  the  classic  VEGF  pathway. 
When  studying  survival  using  the  glioblastoma  cohort  in 
TCGA,  we  demonstrated  a  survival  advantage  correlating  Avith 
the  PDGF-C  expression  only  in  the  classic  subgroup  of  glio¬ 
blastoma.  We  hypothesize  that  the  correlation  Avithin  the 
classic  subgroup  might  be  due  to  differences  of  the  origin  of 
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